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Abstract 

The theory of first strain gradient elasticity (SGE) is widely used to model size and non-local effects 
observed in materials and structures. For a material whose microstructure is centrosymmetric, SGE 
is characterized by a sixth-order elastic tensor in addition to the classical fourth-order elastic tensor. 
Even though the matrix form of the sixth-order elastic tensor is well-known in the isotropic case, its 
complete matrix representations seem to remain unavailable in the anisotropic cases. In the present 
paper, the explicit matrix representations of the sixth-order elastic tensor are derived and given for 
all the 3D anisotropic cases in a compact and well-structured way. These matrix representations 
are necessary to the development and application of SGE for anisotropic materials. 
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1. Introduction 

In classical continuum mechanics (Truesdell and Toupin (1960); Truesdell and Noll (1965)), only 
the first displacement gradient is involved and all the higher order displacement gradients are ne- 
glected in measuring the deformations of a body. This usual kinematical framework turns out not 
to be rich enough to describe a variety of important mechanical and physical phenomena. In par- 
ticular, the size effects and non-local behaviors due to the discrete nature of matter at a sufficiently 
small scale, the presence of microstructural defects or the existence of internal constraints cannot 
be captured by classical continuum mechanics (see, e.g., Marangantia and Sharma (2007) and the 
references cited therein for more details). The early development of high-order (or generalized) con- 
tinuum theories of elasticity was undertaken in the 1960s and marked with the major contributions of 
Toupin (1962); Koiter (1964); Mindlin (1964, 1965); Eringen (1968); Mindlin and Eshel (1968). For 
the last two decades, the development and application of high-order continuum theories have recently 
gained an impetus, owing to a growing interest in modeling and simulating size effects and non local 
behaviors observed in a variety of materials, such as polycrystalline materials, geomaterials, bioma- 
terials and nanostructured materials, and in small size structures (see, e.g.. Fleck and Hutchinson 
(1997); Nix and Gao (1998); Lam et al. (2003); dell'Isola et al. (2009, 2011); Liu et al. (2011)). At 
the same time, the development of homogenization theories (see, e.g., Forest (1998); 
Kouznetsova et al. (2004); Trinh et al. (2012)) makes it possible to determine higher- 
order moduli in terms of material local properties and microstructure. In partic- 
ular, it has been recently evidenced (Alibert et al. (2003); Seppecher et al. (2011)) 
that microstructures can be designed to render higher-grade effects predominant. 
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Prom the numerical point of view, theories of generalized continua, such as strain- 
gradient theory, can be used to avoid explicitly meshing coarse heterogeneous materials 
(Kruch and Forest (1998); Tekoglu and Onck (2008); Pau and Trovalusci (2012)). 

The theory of first strain-gradient elasticity (SGE) proposed by Mindlin (1964) and Mindlin and Eshel 
(1968) is among the most important high-order continuum theories which have been elaborated for 
the last half century, and it is currently widely used. In this theory, the infinitesimal strain tensor 
£ and its gradient u = Ve are linearly related to the second-order Cauchy stress tensor a and the 
third-order hyperstress tensor r by equation (1) where a fourth-order tensor C, a fifth-order tensor 
M and a sixth-order tensor A are involved and verify the index permutation symmetry properties 
(2) and (3). The simplest one of these three tensors, C, defines the conventional elastic properties 
of a material. The study of C had experienced a long history (Love (1944)) before a complete 
understanding of C was achieved quite recently. In fact, only about 20 years ago and for the 
first time, Huo and Del Piero (1991) explicitly posed, rigorously formulated and treated the fun- 
damental problem of determining all the rotational symmetry classes that the fourth-order elastic 
tensor C can possess. This problem has then received the attention of researchers from mechan- 
ics, materials science, physics, applied mathematics and engineering (see, e.g., Zheng and Boehler 
(1994); Forte and Vianello (1996, 1997); He and Zheng (1996); Xiao (1997); Chadwick et al. (2001); 
Bona et al. (2004, 2007); Moakher and Norris (2006)). A comprehensive understanding of C is now 
available in the sense that the correct answers to the following three fundamental questions have 
been provided: 

(a) How many symmetry classes and which symmetry classes has C ? 

(b) For every given symmetry class, how many independent material parameters has C ? 

(c) For each given symmetry class, what is the explicit matrix form of C relative to an orthonormal 
basis ? 

However, these questions remains largely open in regard to the fifth-order tensor M and sixth-order 
tensor A in the theory of SGE. Indeed, in the 3D isotropic situation, Toupin (1962) and Mindlin 
(1965) gave the general form of A containing 5 independent material parameters, and Mindlin (1964) 
showed that M must be zero. In the 2D context, Auffray et al. (2009a) derived all anisotropic ma- 
trices of A. At the present time, in the 3D case, few results are known. A first result was established 
by deirisola et al. (2009), they constructed and studied a matrix representation of the sixth-order 
tensor A in the isotropic situation. And, more recently, Papanicolopulos (2011) investigates fea- 
tures of the M tensor in the S0(3)-symmetry. For that symmetry class a coupling between the the 
gradient and the first gradient exists. But, apart from these results, little is known about M and A. 

In developing and applying the theory of SGE, there is currently a real need for posing and 
answering the foregoing three questions about the fifth-order tensor M and sixth-order tensor A. 
We first consider materials whose microstructure is centro-symmetric. In this case, M becomes zero 
and investigations can be carried out only on A. In a companion paper dedicated to 3D SGE 
Le Quang et al. (2012), we have proved that A has 17 symmetry classes, identified the nature of 
each symmetry class and determined the number of independent material parameters of A belonging 
to a given symmetry class. Nevertheless, in the literature and in our aforementioned work, the 3D 
explicit matrix forms of A have not been furnished for the 16 anisotropic symmetry classes. This 
situation prevents the theory of first SGE from being developed and applied for anisotropic materials. 

Compared with the classical fourth-order tensor C, the sixth-order tensor A is much more 
complex and richer. Indeed, A has 16 anisotropic symmetry classes whereas C possesses 7 ones. 
In this regard, remark that, for example, the transverse hemitropy and transverse isotropy are two 
distinct symmetry classes for A but shrink into one symmetry class for C. A similar phenomenon 
also occurs for the tetrahedral and cubic symmetries. In addition, even for the same symmetry class, 
the number of independent material parameters of A is much higher than that of C. For instance. 
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the cubic symmetry, C contains 3 independent parameters while A comprises 11 ones. 

The present work aims to solve the problem of obtaining the base of the explicit matrix repre- 
sentations of A for all the 17 symmetry classes. As will be seen, the complexity and richness of A 
make a proper solution to this problem not straightforward at all. In fact, even though the matrix 
form of A relative to an orthonormal basis is known, how to express A as a symmetric square matrix 
in a compact and well-structured way is far from being obvious. 

The paper is organized as follows. In the next section, the constitutive law of SGE is recalled 
and the essential results obtained by Le Quang et al. (2012) about the symmetry classes of A are 
recapitulated together with the most important concepts involved. The main results obtained by 
the present work are given in section 3. They include the explicit matrix representations of A 
for the 17 symmetry classes, which are rendered very compact and well-structured by proposing 
an original three-to-one subscript correspondence. Each symmetry class is associated to a simple 
geometric figure, and the matrix representations of A are presented in such a manner that they 
can be directly used without resorting to the theory of rotational groups. In section 4, the logic of 
the three-to-one subscript correspondence chosen in section 3 is explained, the general structure of 
the matrix representations of A is highlighted, and some salient differences between first SGE and 
classical elasticity are emphasized. In section 5, a few concluding remarks are drawn. 



2. Strain-gradient elasticity 

2.1. Constitutive law 

In the (first) strain-gradient theory of linear elasticity (see, e.g., Mindlin (1964); Mindlin and Eshel 
(1968)), the constitutive law gives the symmetric Cauchy stress tensor cr and the hyperstress tensor 
T in terms of the infinitesimal strain tensor e and strain-gradient tensor u) = Ve through the two 
linear relations: 

'Tijk — ^Imijk^lm ~l~ ■^ijklvan^lmn- 

Above, (Tjj, ffij, Tijk and Uijk = Eij^k are, respectively, the matrix components of a, e, t and Ve 
relative to an orthonormal basis {ei,e2,e2} of a three-dimensional (3D) Euclidean space; Cijim, 
Mijimn and Aijkimn are the matrix components of the fourth-, fifth- and sixth-order elastic stiffness 
tensors C, M and A, respectively. These matrix components have the following index permutation 
symmetry properties: 

Cijlm CjiiYTi Cijnij, ^ijklm ^jiklm ^ijlkmi (2) 



-Aijkimn — Ajif^iYnn — AiYnnijk- (3) 

In the case where the microstructure of a material exhibits centro-symmetry, the fifth-order elastic 
stiffness tensor M of this material is null, so that the constitutive law (1) becomes uncoupled: 

In this paper, we are essentially interested in answering the question of what are the possible different 
matrix forms for A with the index symmetry property (3), referred to as the strain-gradient elasticity 
(SGE) tensor. The same question can be posed for M and will be treated in another paper. 
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2.2. Symmetry classes 

In a recent paper (Le Quang et al. (2012)), we have solved the problem of determining all the 
symmetry classes that the sixth-order SGE tensor A can have. For the purpose of the present paper, 
we below recall some relevant definitions and summarize the main result obtained in Le Quang et al. 
(2012). 

Let Q be an element of the 3D rotation group S0(3). An SGE tensor A is said to be invariant 
under the action of Q if 

QioQjpQkqQlrQmsQnt-^opqrst — -^ijklmn- (5) 

The symmetry group of A is defined as the subgroup Ga of SO (3) formed of all the 3D rotation 
tensors leaving A invariant: 

Ga = {Q G S0(3) \ QioQ jpQkqQlrQ msQnt -^opqrst — -^ijklmn }• (6) 

From the physical point of view, it is meaningful to consider two SGE tensors A and IB as exhibiting 
rotational symmetry of the same kind if their symmetry groups are conjugate in the sense that 

there exists a Q G S0(3) such that Gb = QGaQ^- (7) 

Thus, the (rotational) symmetry class associated to an SGE tensor A can be naturally defined as 
the set [Ga] of all the subgroups of S0(3) conjugate to Ga: 

[Ga] = {G C S0(3) |G = QGaQ^, Q gS0(3) }. (8) 

Li other words, the symmetry class to which A belongs corresponds to its symmetry group modulo 
S0(3). In fact, this classification leads to a partition of the set consisting of all the subgroups of 
S0(3). 

For later use, we introduce some additional notations. First, the rotation about a vector a 
through an angle 9 G [0,27r) is denoted by Q(a,0); in particular, the rotations Q(ei + e2 + ea, 27r/3) 
and Q[2(\/5 - l)e2 + (^/5 + l)e3,27r/3] are in short noted as Q for Q, respectively. In what 
follows, use will be made of the subsequent standard group notations: the second-order identity 
tensor group I; the 2D rotation group S0(2) consisting of all rotations Q about a 3D vector, say 
e3, such that Qes = 63; the 2D orthogonal group 0(2) composed of all orthogonal tensors Q 
such that Qe3 = ±63 for a fixed direction 63; the cyclic group with r > 2 elements generated 
by Q(e3,27r/r); the dihedral group D^ (r > 2) with 2r elements generated by Q(e3,27r/r) and 
Q(ei,7r); the tetrahedral group T with 12 elements generated by D2 and Q; the octahedral group 
O containing 24 elements generated by D4 and Q; the icosahedral group X having 60 elements 
generated by D5 and Q. Recall that the subgroups T, O and X map a tetrahedron, a cube and an 
icosahedron onto themselves, respectively. 

In the recent paper of Le Quang et al. (2012), it is proved that the number of all possible sym- 
metry classes for the SGE tensors is 17. In addition, the number of independent matrix components 
of an SGE tensor belonging to a given symmetry class has also been determined by Le Quang et al. 
(2012). These results are summarized in Table 1 for the purpose of the present work. And all the 
17 symmetry classes are graphically illustrated in Figures 1 to 17. 

Note that the number of all possible symmetry classes for the SGE tensors is much higher 
than the one for the classical elasticity tensors, since the former is 17 while the latter is 8. In the 
totally anisotropic case where [Ga] = I and in the isotropic case where [Ga] = SO (3), the number of 
independent components of A is equal to 171 and 5, respectively. This clearly shows the SGE theory 
is much more complicated than the classical elasticity theory where the corresponding numbers are 
21 and 2, respectively. 
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Name 


Triclinic 


Monoclinie 


Orthotropic 


Chirally trigonal 


Trigonal 


[Ga] 


I 


[Z2] 


[D2] 


[Z3] 


[D3] 


#mdcp(A) 


171 


91 


51 


57 


34 


Name 


Chirally tetragonal 


Tetragonal 


Chirally pentagonal 


Pentagonal 


Chirally hexagonal 


[Ga] 


[Z4] 


[D4] 


[Z5] 


[Ds] 


[Ze] 




45 


28 


35 


23 


33 


Name 


Hexagonal 


Transversely hemitropic 


Transversely isotropic 


Tetrahedral 


Cubic 


[Ga] 


m 


[S0(2)] 


[0(2)] 


m 


[O] 


#mdcp(''^) 


22 


31 


21 


17 


11 


Name 


leosahedral 


Isotropic 








[Ga] 




S0(3) 








#indcp(A) 


6 


5 









Table 1: The names, the sets of subgroups [Ga] and the numbers of independent components #indcp(A) for the 17 
symmetry classes of SGE. 



3. Matrix representations of strain-gradient elasticity: main results 

For an SGE tensor A belonging to one of the 17 symmetry classes listed in Table 1, it is very 
important in various situations of theoretical and practical interest to know the explicit matrix form 
of A relative to an orthonormal basis {61,62,63}. To solve this problem in a satisfactory way, we 
need: (i) firstly to identify what are the non-zero matrix components of A and what are the possible 
relations between its non-zero components; (ii) secondly to elaborate an appropriate representation 
method according to which the matrix of A is well-structured and takes a simple and compact form 
so as to be used easily. The solution to the first part of the problem can be found in the paper of 
Le Quang et al. (2012). The solution to the second part of the problem is still lacking and will be 
provided in what follows. 

3.1. Orthonormal basis and matrix component ordering 

Let be defined the following subspace of third-order tensors 

3 

5^ = {T|T= Tykei0ej0ek, T,jk = T^ik} (9) 

i,j,k=l 

which is an 18-dimensional vector space. By (4), an SGE tensor A is a linear symmetric transfor- 
mation from to S^. In order to express the strain gradient a; or the hyperstress tensor r as 
a 18-dimensional vector and write A as a 18 x 18 symmetric matrix, we introduce the following 
orthonormal basis vectors: 



V2 



+ Y^(^i^^j+^j'^^i)'^^k, 1 < a < 18 (10) 



where the summation convention for a repeated subscript does not apply. Then, the strain gradient 
u, the hyperstress tensor r and the SGT A can be expressed as 

18 18 18 

a; = ^0)060,, T = ^fQ,ea, A= ^ AapGa'^ep, (11) 

a=l a=l a, 13=1 
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-i 




o 
o 


/I 







ijk 


1 1 -i 
111 


221 


122 


331 


loo 

133 


Frivileged direction: 1 


a 


6 


7 


8 


9 


10 




ijk 


222 


112 


121 


332 


233 


Privileged direction: 2 


a 


11 


12 


13 


14 


15 




ijk 


333 


113 


131 


223 


232 


Privileged direction: 3 


a 


16 


17 


18 








ijk 


123 


132 


231 






No privileged direction 



Table 2: The three-to-one subscript correspondence for 2D (in boldface in the table)and 3D strain-gradient elasticity 



SO that the second SGE relation in (4) can be conveniently written in the matrix form 

Ta = AapOjp. (12) 

Note that, with the orthonormal basis (10), the relationship between the matrix components uja 
and ojijki the one between fa and Tj^-fc, and the one between A^p and Aij^imn are specified by 

co^ = hl^ '='^ fa = i^f'''='^ (13) 

IV2a;ijfc if z/j; [v^Tijk if z / j; 



An 



Ajklmn if i=j and I = m, 

V^-^ijklmn if ^ 7^ J and I = m or i = j and / 7^ m, (14) 
2Aijfc/mn if j and / / m. 



It remains to choose an appropriate three-to-one subscript correspondence between ijk and a. 
For the SGE matrix A^is to be well-structured and exhibit a simple compact form for a given symme- 
try group Ga, some criteria guiding the choice of an efficient three-to-one subscript correspondence 
are necessary. The criteria adopted in this work are detailed and explained in the next section. The 
final three-to-one subscript correspondence is specified in Table 2. 

3.2. Rotation matrix 

With the basis defined in (10) and the three-to-one subscript correspondence detailed in Table 
2, the action of a rotation tensor Q G S0(3) on an SGE tensor A can be represented by a 18 x 18 
rotational matrix Q in such a way that 

QioQjpQkqQlrQmsQnt-^opqrst — 

where ^ 

Qa/S ~ '^i^QioQ jp ~\~ QipQ jo)Qkq (16) 

with Q and /3 being associated to ijk and opg, respectively. Thus, formula (5) expressing the 
invariance of A under the action of Q is equivalent to 

QAQ^ = A (17) 

where Q stands for the 18 x 18 matrix of components Qap and A the 18 x 18 matrix of components 
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Figure 1; Triclinic system (I-invariance): the material is completely asymmetric. 



3.3. Matrix representations for all symmetry classes 

We are now ready to give the explicit expression of the SGE matrix A for each of the 17 symmetry 
classes. To be expressed in a simple and compact way, the matrix A for every symmetry class is 
split into sub-matrices so as to make appear elementary building blocks. The order adopted to 
specify the expressions of A for the symmetry classes [Z^] and [D^,.] is A; = 1, 2, 4, 3, 6, 5 and oo. 
The reason for such a choice is explained in subsection 4.2. 

3.3.1. Symmetry class characterized by I 

In this most general case, illustrated by figure 1, the material in question is totally anisotropic 
and the SGE matrix A comprises 171 independent components. The explicit expression of A as a 
full 18 X 18 symmetric matrix can be directly obtained by formula (14). We first define 



the 



n(n+l) 



dimensional space Ai^{n) consisting of n x n symmetric matrices; 



• the n^-dimensional space Ai{n) made of n x n matrices; 

• the nm-dimensional space M{n,m) composed of n x m matrices. 

Then, we can write A in the following way 



£;(15) f{25) q(15) 
/f(15) /(15) 



where the subscript S indicates that the matrix is symmetric and the form and number of indepen- 
dent components of each involved sub-matrix are specified by 

. i?(15), e A^5(5). 

. 5(25), ^(25), f(25) e _^(5). 

• ^(15), /(15) g7W(5,3); 

. J(6) G 7W^(3). 

For example, ^^^^^ is an element of 7W'^(5) and contains 15 independent components while S^^s) 
belongs to A1(5) and comprises 25 independent components. 
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A ^2 




Figure 2: Monoclinic system (Z2-invariance); the material is vr- invariant about eg. 

Ae2 





5 ^ 


63" 


ei 



Figure 3: Orthotropic system (D2-invariance): the material is 7r-invariant about each of ei, 62 and 63. 



3.3.2. Symmetry classes [Z2] and [D2] 

According as the symmetry class [Z2] or [D2] is concerned (c.f. figures 2 and 3), the mate- 
rial described is monoclinic or orthotropic, and the SGE matrix A contains 91 or 51 independent 
components. Using the three-to-one subscript correspondence given in Table 2, the SGE matrices 
exhibiting the Z2-symmetry and Z)2-symmetry are well-structured and take the compact forms: 



^^(15) ^(25) 

H(15) /(15) 



V 





^(15) p(25) Q 
J/(15) 



A 



D2 





£(15) 







iy(i5) 



\ 



where ^(1^), ^(15)^ ^(15) g M^{5), 5(25),f(25) e 7^(5), /(is), e A^(5, 3) and e A^^(3). 

Here, because of its practical interest we give two conjugate representations of the same sym- 
metry class [Z2]. In the first representation the vr-rotation is taken around es as indicated by the 
notation Z2^, meanwhile in the second case the rotation is taken around ei (as indicated by Z2^). 
The first situation is considered in order to be coherent with the representation of the other cyclic 
classes, in which the generating rotation is taken around 63. The second representation we exhibit 
correspond to the common case of a monoclinic material, the combination of the Zg^-invariance and 
the central inversion (always contained in the symmetry group of any even-order tensor) leads to 
the existence of symmetry plane which normal is 63. 

It is remarkable that the non-zero matrix blocks of Aj)^ are diagonally located. Note that A^b^ 

and A^'^s are identical to Ad^ to within the two non-diagonal matrix blocks B^^^) and in one 
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Figure 4: Chirally tetragonal system (Z4-invariance): the material is ^--invariant about 63. 





'62 














► 

ei 



Figure 5: Tetragonal system (D4-invariance) : the material is ^-invariant about ea, and 7r-invariant about 
(Q(e3,fc7r/4)ei) with fc G Z. 

case and to within £'(2^) and D^^^^ in the other. 

3.3.3. Symmetry classes [Z4] and [D4] 

The materials characterized by the symmetry classes [D4] and [Z4], shown figures 4 and 5, are 
said to be tetragonal and chirally tetragonal. The numbers of independent components of A with 
the Z)4-symmetry and Z4-symmetry are 28 and 45, respectively. To write the corresponding SGE 
matrices in a compact way, we first introduce 

• the "^"2 ""^^ -dimensional space A^^(n) consisting of n x n anti-symmetric matrices; 

• the matrices H^^\ /^^^ and J*^^^ with 9, 7 and 4 independent components defined by 



/ hii hi2 hi3 hi2 hi3\ 

h22 h23 h2A h25 

h33 h25 /l35 

^22 ^23 



1(7) 



( ii2 -«12\ 

«21 «22 «23 

«31 *32 *33 

-«21 — i23 — *22 

\-«31 -433 — «32 



J(4) 



J12 J12 

i22 j23 
i22 y 
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Figure 6: Chirally trigonal system (Zs-invariance): the material is ^-invariant about 63. 




Figure 7: Trigonal system (Ds-invariance) 
(Q(e3,fc7r/3)ei), A; € Z. 



the material is 



2tt_ 
3 ■ 



invariant about 63, and vr-invariant about 



Then, the Z4-symmetric and L'4-symmetric SGE matrices can be written as 



Z4 



SdO) \ 

A(i^) 
_H-(9) 7(7) 



^D4 
























V 







where A^^^) G {f>) and S^^") e A^"^(5). Owing to the subscript ordering specified by Table 2, 
both the expressions of Az_^ and ^4^4 exhibit compact structure. In addition, Ay)^ has a diagonal 
block structure. 

3.3.4- Symmetry classes [Z3] and [D3] 

The materials having the symmetry classes [D3] and [Z^], shown figures 6 and 7, are referred 
to as being trigonal and chirally trigonal, respectively. The numbers of independent components 
contained in the corresponding matrices Az-^ and Ad^ are 57 and 34, respectively. As will be seen, 
even if use is made of the subscript ordering of Table 2, the matrix expressions of Az^ and ^d3 
remain quite complex. However, it is possible to get a good understanding of the structures of 
Az^ and Ad^ by defining appropriate independent sub-matrix blocks and making appear dependent 
sub-matrix blocks. More precisely, Az^ and Ad^ can be expressed as 



A' 



Z3 



G(9) 



+ 



/(G(9)) 

/(Z)(4)) 



f(FW)\ 



/ 



10 



A 



D3 



+ 



/o 

/(DW) 

/(JW) 

V / 



/a(ii)+7?Ac DWNy 

^(11) 


In these two expressions, rj and are two scalar material parameters; A^^^\ B^^\ C^'^\ D^^\ F^^\ 
G(9),ii'(6),/(4) and J(4) are 9 independent sub-matrices; Ac and Be are two coupling matrices con- 
taining no material parameters; /(G^^^), /(F^^)), /(Z?(^)) and /(J(^)) are the matrix-value functions 
of G(9), and J^^\ respectively. 

First, the expression of A^^^^ with 11 independent components is specified by 



(a 



^(11) 



11 112 
0122 



ai3 ai4 ais"^ 

-ai3+V2oni "I "11 

— ai2+ani '^34 ass 

a44 a45 

0-55 



where 



ai = ai4 



\/2( 



034 , «ii = ai5 



035 5 "III 



an - 022 



«iii 



Next, the expressions 

of 5(6), C(3), Z)(4), and J^^) are given by 

/o 61 



012 




'^bl2 624 + V2634 b25 + V2635\ 

_V2 







^24 
634 





v 



fe25 
635 
645 





an + 022 



C(3) 



/ cii 

— Cll 

-v^cii 





i7(8) 



with 



Cl2 
-Cl2 






/ hi 

-hi 

-V2hi 




V 



Cl3 
-Cl3 




/12 

-/l2+/3l 
-V^{/l2 






C12 

-Cl2 






Cl3 \ 

— Cl3 
-V2C-13 





/l3 

/23 

-v^(/i5+/3ii) 



Z)(4) 



/14 

-/l2+/3l 
-V^{/l2-4 





\d5i 



Pi 



^(9) 





/43 







/ss 





2 - fu 

2 




_ /l3 + /23 

2 


911 


912 


913 


921 


923-27ni 


923 


V271 


V2711 


V2(27iii+7ii) 


941 


942 


942 


951 


952 


952 



<ii2 
-di2 
-V2di2 



/l5 

-/i5-2/3ii 

^V^(/l3-/9ll) 
-/43 
— /53 



-rfl2\ 

dl2 
V2di2 






/ 



11 



with 



^(6) 



7i -- 

All 

V 



511 - 521 



hi2 
h22 



2 




111 






hi2 


hi3\ 




h23 


h22 


h23 




h33 


h23 


h33 






h22 


h23 








h33/ 


S 



913 - 923 



7iii 



512 - 513 



j(4) 



/ 





ii2 


-ii2 \ 







«22 


— «22 — ^2*31 




«31 


«32 


—i32 


V 





122 + ^2131 


-i22 


-hi 


«32 


-«32 / 



The matrices Ac and Be are independent of material parameters and take the following forms: 
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Finally, the matrix- value functions f{G^^^), f{F^^')), f{D^^^) and /(J^^-*) are defined by 
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Figure 8: Chirally hexagonal system (Zo-invariance): the material is -j-invariant about ea 




Figure 9: Hexagonal system (De-invariance); the material is ^-invariant about ea, and vr-invariant about 
(Q(e3,fc7r/6)ei), k € Z. 



A^2 




Figure 10: Chirally pentagonal system (Zs-invariance): the material is ^-invariant about 63 



Ae2 




Figure 11: Pentagonal system (Ds-invariance): the material is ^-invariant about es, and 7r-invariant about 
(Q(e3,fc7r/5)ei), fc e Z. 



3.3.5. Symmetry classes [Zg] and [Dq] 

The hexagonal and chirally hexagonal materials are described by the symmetry classes [Dq] and 
[Zq], illustrated by figures 8 and 9, and have 22 and 33 independent parameters, respectively. They 
can be considered as being degenerated from the trigonal and chirally trigonal materials. Precisely, 
the associated SGE matrices Azg and Aj)^ take the following simpler forms: 
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where rj and 6 are two scalar parameters and the non-zero matrix blocks are identical to the relevant 
ones given for the symmetry class [Z3]. 

3.3.6. Symmetry classes [Z^] and [D^]. 

According as a material belongs to the pentagonal symmetry class [D5] or chirally pentagonal 
symmetry [Z^], its number of independent parameters is 23 or 35. The corresponding SGE matrices 
(see figures 10 and 11) Az^ and Ay)^ are given by 
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In these two expressions, the sub-matrices A^^^\ B^^\ H^^\ j(4) 

are specified in the 

foregoing case where the trigonal and chirahy trigonal symmetry classes, [D3] and [Z3], are concerned. 
The remaining sub-matrices F^^^ and G^^\ each of which contains 2 independent components, and 
the matrix-value functions /(F^^^), /(G^^^) and g(J^^^) take the followings forms: 
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5.5.7. Symmetry classes [S0(2)] and [0(2)]. 

These two symmetry classes, shown by figures 12 and 13, characterize the transversely hemitropic 
and transversely isotropic materials, respectively. The associated SGE matrices j4go(2) and ^0(2)1 
containing 31 and 21 independent components, respectively, have the following expressions: 
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H(6) /(4) + J(J(4)) 

\ /o 



~^ /(JW) 

where the sub-matrices A^^^\ B^^\ H^^\ I^^^ and j(^) are defined in studying the trigonal and 
chirally trigonal symmetry classes [D3] and [Z3]. 
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Figure 12; Transversely hemitropic system (S0(2)-invariance): the material is oo-invariant about es 





Figure 13: Transversely isotropic system (0(2)-invariance): the material is oo-invariant about 63, and 7r-invariant 
about any in-plane axis. 
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ei 

V 



Figure 14: Tetrahedral system (T-invariance): the material is D2-invariant about ea and invariant under the permu- 
tation (ei — 5- 62 — > e-j — > ei) (The tetrahedron is drawn inside the cube represented in fig. 15). 



3.3.8. Symmetry classes [T] and [O]. 

The materials described by the symmetry classes [T] and [O] shown by figures 14 and 15, 
are said to be tetrahedral and cubic, respectively. The former is characterized by 17 independent 
material parameters and the latter by 11 ones. The corresponding SGE matrices Aj- and Aq have 
the expressions: 
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where A^^^^ is an element of M^{5), P is the permutation matrix defined by 
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and the sub-matrices A^^^ and J'-^^ with 9 and 2 independent material parameters, are specified by 
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Note that both Aj- and Aq have a very compact diagonal structure. 
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Figure 15: Cubic system (O-invariance): the material is D4-invariant w.r.t 63, and invariant under the permutation 
(ei — )• 62 — >■ ea — >■ ei). 




Figure 16: Icosahedral (I-invariance): the material is T-invariant about 63, and Zs-invariant about v = (e2 + (l — 0)63). 
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3.3.9. Symmetry classes [I] and [S0(3)] 

The last symmetry classes [I] and [S0(3)] shown by figures 15 and 16, are icosahedral and 
isotropic. The corresponding SGE matrices Ax and ^sO(3) comprises 6 and 5 independent material 
parameters, respectively. They have the following expressions: 
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where r/ is a scalar material parameter, the 5x5 sub-matrix J^^^ contains 5 independent material 
parameters, P is the same permutation matrix as the one defined in treating the tetrahedral sym- 
metry class [7^, is a 5 X 5 sub-matrix containing no material parameter, Jc is a 3 x 5 sub-matrix 
involving the the golden number and /(^^^^)is a matrix- value function of J^^^ . 
More precisely, AS^^ takes the form 



with 



A>j^ and Jc are given by 
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with (f) being the conjugate of the golden number (j) defined as = = ! — (/), and f{A^^^) is 

specified by 

«V+V2aiv aill-a^v 

with 

av = «22 — 012- 



4. Matrix representations of strain-gradient elasticity: some general remarks 

In the previous section, we have presented the main results about the matrix representations 
for the 17 symmetry classes of strain- gradient elasticity. These matrix representations have a very 
compact structure and exhibit some general properties. In this section, we explain the reasons 
underlying the three-to-one subscript correspondence specified in Table 2 and make some general 
remarks. 



Matrix component ordering 

In subsection 3.1, starting from a 3-dimensional orthogonal basis {ei,e2,e3}, we have con- 
structed an 18-dimensional orthonormal basis {ei,e2, ...,ei8} for S^. An SGE tensor A is a sym- 
metric linear transformation from to . For its matrix representation relative to the basis 
{ei,e2, ...,ei8} to be well-structured, some criteria have to be established to make a good choice of 
the three-to-one subscript correspondence between ijk and a. The criteria we have elaborated are 
explained below. 

First, we consider a cubic material which is characterized by the octahedral group O graphically 
illustrated by figure 15. In this case, the three-to-one subscript correspondence between ijk and a 
is required to be such that: 

(i) the matrix Aq is block-diagonal; 

(ii) each diagonal block matrix of Aq contains no zero components; 

(iii) each diagonal block matrix of Aq is invariant under every cyclic permutation of ei, 62 and 
63- 

Next, we are interested in a tetragonal material described by the tetragonal group D4 described by 
figure 5. In this situation, we require the three-to-one subscript correspondence to be such that 

(iv) the diagonal block matrices of A-q^ related to the plane ei — e2 are invariant under the 
permutation of ei and e2. 

The satisfaction of the foregoing requirements (i)-(iv) has the consequence that Aq and A-q^ 
take the forms 
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with A € A4^{5), H G A4^(5) and J G A4^{3). In fact, the requirement (i) gives the general 
shape of Ao but does not fix the number of block matrices of Aq. The first idea, which seems 
"natural", is to decompose the diagonal part of Aq into three 6x6 block matrices. However, use of 
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this decomposition makes appear some zero components in each block matrix. The elimination of 
zero components inside each diagonal block matrix motivates the requirement (ii) and is performed 
by carrying out column/row permutations, leading to the decomposition the diagonal part of Aq 
or into three 5x5 block matrices plus one 3x3 matrix. The requirements (iii) and (iv) are 
destined to order the columns and rows within each block matrix. Precisely, the condition (iii) leads 
to the invariance of every block matrix under a cubic symmetry transformation. In particular, the 
first three diagonal matrices are identical. The condition (iv) is imposed for the first two diagonal 
block matrices of to be the same and for the privileged axis defined by to be distinguished 
from the privileged axes defined by ei and 62 which share the same symmetry status. 

The three-to-one subscript correspondence between ijk and a specified by Table 2 is established 
in agreement with the foregoing requirements (i)-(iv). In Table 2, the second, fourth and sixth 
rows are schemed out by singling out one privileged direction and describe the interactions of the 
remaining directions with the privileged one. The eighth row is "mixed" in the sense that it involves 
all the three directions. 

Note that the fourth row is deduced from the second row by the transposition (12) while the sixth 
row is obtained from the second row by an anti-cyclic permutation (132). The reason for doing so 
instead of deducing the fourth and sixth rows from the second row by an cyclic permutation (123) 
is twofold. First, since most of the symmetry classes of SGE are planar, it appears judicious to 
privilege the planar symmetry classes, i.e., S0(2), [Z^] and [D^] with A: = 2, 3, 4, 5 and 6, in 
structuring the matrix representations of SGE. Next, our 3D matrix representations of SGE can be 
easily degenerated into the 2D ones by conserving only the elements of A^is with a,/3 = 1, 2, 3, 6, 
7 and 8 (see the up left elements in Table 2). 



4-2. Generic matrix forms 

As a consequence of the foregoing three-to-one subscript correspondence, the matrix represen- 
tations of SGE have the following generic forms for the different symmetry classes. More precisely, 
for the chiral planar symmetry classes [Z^] with A; = 2, 3, 4, 5 and 6, the matrix representations 
exhibit the generic shapes 



/ Di E12 Ei3 Cl^ 
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\ J J 
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V J J 



s \ " / s 

Concerning the dihedral symmetry classes [D^] with /c = 2, 3, 4, 5 and 6, we have 
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\ JJ 



The matrix representations for the spatial symmetry classes [T], [O], [I] and S0(3) take the generic 
form 

/DO 0\ 

V{D) 
D 

V J/. 

where V means a permutation of the block matrix D. In particular, for the symmetry classes [O] 
and S0(3), the permutation reduces to identity, so that V{D) = D. 



As 
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From the above generic matrix forms for the different symmetry classes, it can be seen that the 
elementary blocks involved in the 3D matrix representations of SGE are of four types: 

(a) D-type block diagonal matrices belonging to A^'^(5); 

(b) J- type block diagonal matrices belonging to A^'^(3); 

(c) E'-type block extra-diagonal matrices belonging to A4{5); 

(d) C-type block extra-diagonal matrices belonging to A^(5,3). 

4-3. Remarkable differences between SGE and classical elasticity 

As recalled in section 2, SGE has 17 symmetry classes while classical elasticity possesses 8 
symmetry classes. 

• First, note that the non-crystallographic symmetry classes [Z5], [D5], [T] and [I] make sense in 
the case of SGE but disappear in classical elasticity. These new classes are meaningful for the study 
of quasi-crystallographic alloys. As is well-known, most of quasi-crystallographic ordered materials 
exhibit icosahedral symmetry (Gratias et al., 2000). In addition [Z5] and [D5] symmetry classes 
are related to Penrose tilling (Penrose, 1974), a well- know toy- model to understand properties of 
quasi-crystallographic materials . 

• Next, SGE is sensitive to chirality: (i) the chiral symmetry classes [Z3] and [Z4] are present in 
SGE but absent in classical elasticity; (ii) the transversely hemitropic symmetry class [SO (2)] holds 
for SGE but not for classical elasticity. 

• Last, the hexagonal and chirally hexagonal symmetry classes [Dg] and [Zg] are meaningful in 
the case of SGE but have the same effects as the ones of [0(2)] in the case of classical elasticity. 
The combination of this result with the crystallographic restriction theorem^ leads to the fact that 
a 2D periodic medium can neither be transversely isotropic nor transversely hemitropic. To sum 
up, both for the 2D and 3D cases, every periodic material is anisotropic for SGE. 

The aforementioned differences constitute one of the reasons for which the 3D matrix represen- 
tations of SGE are much more complex and subtler than those of classical elasticity. 

4.4- The chiral sensitivity of SGE: comments and numerical example 

As pointed out above, SGE is sensitive to chirality. Indeed, there exist two different types of 
chirality for SGE, which are explained below. 

• The first type of chirality is related to the chiral subgroups of 0(2), i.e., SO(2)-subgroups. 
This is an "in-plane" chirality which couples spatial directions, and will be called 5-type. In 
addition, the chirality of 5-type can be encoded both by even- and odd-order tensors. The 
chirality of the sixth-order SGE tensor A studied in the present paper is of 5-type. 

• The second type of chirality is related to the chiral subgroups of 0(3) (for a detailed discussion 
on these subgroups, see, for example, Sternberg (1995)), and will be qualified as O-type. In 
contrast with the chirality of 5-type, the one of O-type couples the first- and second-order 
effects, i.e. the stress depends on the strain-gradient and the hyper-stress on the strain. More- 
over, it can be encoded only by odd-order tensors. The chiral effect of the 5th-order tensor M 
involved in the constitutive law (1) of SGE has been studied, for example, by Papanicolopulos 
(2011) considering a particular case. It should be noted that for some symmetry classes both 
phenomena appear. In such a case the couplings are of SO-type. 



In the 2D case, the crystallographic restriction theorem states that the orders of rotational invariance compatible 
with periodicity are restricted to 1, 2, 3, 4 and 6. 
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To illustrate some physical implications of the chirality of 5-type, an example of homogenization 
incorporating strain gradient effects is now numerically studied. We consider three 2D unit cells: (i) 
the first one shown by Fig. 18 is orthotropic and belongs to the symmetry class [D4]; (ii) the other 
two ones described by Fig. 19 and Fig. 20 are chirally orthotropic and falls into the symmetry class 
[Z4]. To distinguish the unit cell of Fig. 19 from that of Fig. 20, the former is said to be "levogyre" 
while the latter is qualified as being "dextrogyre". Observe that they can be obtained from each 
other by the reflection with respect to the horizontal or vertical middle line. Each cell is a square 
of length 2 m centered at the origin, the slots are rectangles of height 0.2 m and width 0.6 m. In 
the orthotropic case the the upper slot is centered at (0,0.7), and the other ones are obtained by 
rotation of ^. For the levogyre cell, the center of the upper slot is shifted by -0.2 m along xi, 
meanwhile for the dextrogyre the shift is of 0.2 m along the same vector. The remaining slots are 
then obtained by symmetry operations. 




Figure 18: An orthotropic cell 




Figure 19: A chirally orthotropic levogyre cell 




Figure 20: A chirally orthotropic dextrogyre cell 

The choice of a 2D example instead of a 3D one is due to the fact that a 2D example is much 
simpler from the numerical standpoint but suffices for illustrating the sensitivity of SGE to the 
chirality of 5-type in a striking way. In the 2D case, the 3-to-l correspondence to be used is that 
given by the bold characters in the left up part of Table 2. 



23 



To compute the components of the sixth-order SGE tensor A, we apply the finite element 
method (FEM) to the aforementioned three unit cells and prescribe quadratic boundary conditions 
(QBCs)^. Such boundary conditions were proposed in Gologanu et al. (1997) and Forest (1998) 
and further discussed by Auffray et al. (2010). The determination of the components of A consists 
in computing the area averages of the first moments of the stress field in a unit cell produced by 
elementary QBCs. This is an extension of the classical computational homogenization procedure 
for elastic heterogeneous media. 

Let us first consider the orthotropic cell with 4 rectangular slots (Fig. 18). The material forming 
the solid part is a linearly elastic isotropic material with Young's modulus E = 200 GPa and 
Poisson's ratio = 0.3. Its relative density, i.e. the ratio of the area of the porous part to the 
area of the solid part, is p = 0.84. Applying the FEM and the homogenization procedure with 
appropriate QBCs, we obtain^ 
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This matrix is block-diagonal and the two diagonal blocks in it are equal to each other. Such a 
matrix shape is in agreement with the results presented in §3.3.3. 

Next, we consider the chirally orthotropic cells of Fig. 19 and Fig. 20. In these cells, the material 
forming the solid part is identical to the one for the cell of Fig. 18 but the position of each slot is 
changed. Using the same computational method as before, we obtain 
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In contrast with the matrix ^4^4, the matrices A-^l and A^d are no more block-diagonal since 
an antiymmetric block matrix coupling the directions 1 and 2 occurs. Further, even though the 



^As indicated in Forest and Trinh (2011), QBCs leads to an overestimation of A. But, as shown in Auffray et al. 
(2010), the corresponding results are qualitatively correct. 

^The unit of the components is MPa.mm^. Their values were calculated by taking into account the relative density 
of the unit and rounded up. 



24 



diagonal blocks of A^l are identical to those of A^d , the out-of-diagonal antisymmetric blocks of 
A'^L are however different from those of A^^d by a sign. 

The comparison between the matrices ^045 A^l and A^d clearly shows the chiral sensitivity 
of SGE. To get more insight, let us examine the first column of A in detail. The elements of this 
column are determined with the aid of the following QBC 

\Ul{xi,X2) = \x\ 
\u2{xx,X2) = 

imposed on the boundary dVl of a unit cell $7 and by computing the induced hyperstress components 
via ^ 

Tijk = -< {(TijXk + (TikXj) > . 

In this expression < • > is the average operator defined by 

where denotes the apparent area of O and the relative density J) allows the correction of the 
usual average operator due to the presence of the voids in a unit cell (Zybell et al. (2009)). The 
components Tm, T122 and T212 are all non-zero for any of the three cells in question. At the same 
time, the components T211 and T121 are not null only when a chiral unit is concerned. To see the 
last point, we write the explicit expressions allowing the computation of T211 and T121: 
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1 = ^ < (0-12x1 + 0-12x1) >=< auxi >, ri2i = ^ < (0-12x1 + 0-11x2) > 



where the moments of 0-12X1, and 0-112:2 are involved. By examining the relevant fields obtained 
through FEM, it is seen that: (i) for the D4-invariant cell, the field CJ12 is symmetric with respect 
to the middle vertical line (Fig. 21), so that 0-12X1 is well equilibrated (Fig. 22) and the area average 
of 0-12X1 is null; (ii) for a chiral cell, for example, an orthotropic levogyre cell, the field CJ12 has not 
the symmetry relative to the middle vertical line (Fig. 23) and the area average of the field 0-12X1 is 
no more null. In this sense, we can say that the non-zero chiral components of A come from the 
lack of a reflection symmetry of the chiral unit cell in question. 



5. Concluding remarks 

Up to now, the development and application of (first) strain gradient elasticity (SGE) have been 
almost exclusively confined to the isotropic case. The complexity and richness of anisotropic SGE 
remain scarcely exploited. In the present work and the companion one (Le Quang et al. (2012)), we 
have studied materials whose microstructure exhibits centrosymmetry. For these materials, SGE is 
defined by a sixth-order elastic tensor A in addition to the conventional fourth-order elastic tensor C. 
In the companion work (Le Quang et al. (2012)), the tensor A has been shown to have 16 anisotropic 
symmetry classes apart from the isotropic symmetry class. In the present work, the explicit matrix 
representations of A have been presented for all the anisotropic symmetry classes and written in a 
compact and well-structured way. These results will be with no doubt useful for the experimental 
identification, theoretical investigation and numerical implementation of anisotropic SGE. 

In the general case where the microstructure of a material does not exhibit centrosymmetry, a 
fifth-order tensor M intervenes in addition to C and A. Apart from a paper of Papanicolopulos 
(2011) in which the author studied M in the S0(3) symmetry class, the questions concerning the 
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oXY (Achiral) 




L 1 3949X10^'^ 



I 



-1.2 -1 -O.a -0.6 -0.4 -0.2 0.2 0.4 0.6 O.B 1 1.2 ▼ -1.3952x10' 



Figure 21: The stress field ai2 in the cell which is not chiral. 



oXY*X (Achiral) 




A 5.6165x10" 
X10"= 



-1 -0.5 



▼ -5.6253x10" 



Figure 22: The stress moment field ai2Xi in the cell which is not chiral. The area average of this field is null. 
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pXY'X (Levogyre) 




A 7 7698x10" 



' -5.1614x10' 



Figure 24: The stess moment field ai2Xi in the levogyre chiral cell. The area average of this field is not null. 
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symmetry classes and complete matrix representations of M are still entirely open. A natural 
continuation of the present work and the companion one (Le Quang et al. (2012)) will consist in 
finding appropriates answers to these questions relative to M. 

Finally, we remark that, even though the theory of SGE was proposed about a half century ago, 
its development and applications in anisotropic cases are still at their beginning. This is particularly 
pronounced when the fifth-order elastic tensor M is involved. 
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